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摘 要 : 掌握 天 山北 坡 土 壤 水 刨 分 异 的 空间 规律 及 其 驱动 力 , 对 生态 预警 和 土壤 侵蚀 防治 具有 重要 意义 ,为 天 山 


北 坡 区 域 生 态 环境 的 综合 治理 提供 理论 依据 和 数据 支持 。 以 天 山北 坡 中 段 山区 为 例 , 基 于 中 国土 壤 流失 方程 (Chi- 


nese Soil Loss Equation ,CSLE) , 采 


] 野 外 调研 .地 理 信息 系统 数理 统计 和 地 理 探测 器 等 方法 ,定量 分 析 2000 一 2018 


年 研究 区 土壤 水 力 侵蚀 的 时 空格 局 特征 (面积 ,强度 和 地 理 分 布 ) ,借助 地 理 探测 需 探 究 降雨 .地形 .土壤 .植被 对 土 


壤 水 力 侵蚀 强度 的 内 在 驱动 力 。 


结果 表明 :(1) 2000 一 2018 年 天 山北 坡 中 段 山 区 土壤 水 力 侵蚀 强度 主要 以 微 度 、 轻 


度 侵蚀 为 主 ,分 别 占 总 面积 的 32.34%-~40.87% .33.36%~43.01%。 近 20a 的 微 度 . 轻 度 侵蚀 的 面积 均 呈 下 降 趋 势 
(-26.70 km2.a -77.47 km°- a`) ,而 其 他 侵蚀 强度 的 面积 均 呈 上 升 趋势 (22.10~30.96 kim? a!) ,总 体 上 土壤 水 力 侵 蚀 


强度 处 于 增加 趋势 .(2) 整体 土壤 侵蚀 模 数 呈 乌 鲁 木 齐 市 > 


河 子 市 。 天 山北 坡 中 段 山 区 侵蚀 强度 的 空间 分 布 与 降雨 .地 形 土壤、 植被 密切 相关 ,土壤 类 型 为 棕 钙 土 EE, 
BEATE + ,植被 覆盖 度 小 于 15% 坡度 大 于 15° 和 降雨 量 在 400~450 mm 范围 内 的 地 区 为 高 风险 侵蚀 区 域 。(3) 分 异性 


昌吉 市 > 旱 康 市 > 呼 图 壁 县 > 玛 纳 斯 县 > 沙 湾 县 > 石 


的 大 小 由 因子 探测 如 中 g 值 来 度量 ,其 y 值 越 大 ,该 影响 因素 对 土壤 侵蚀 空间 分 布 的 解释 力 越 强 , 旦 降雨 (0.49) > E 
二 类 型 (0.17) > 坡度 (0.11) > 植被 覆盖 度 (0.10)。 不 同 影响 因素 通过 交互 作用 明显 增强 土壤 侵蚀 空间 异 质 性 ,及 植 


被 覆盖 度 与 降雨 因子 的 耦合 作用 ,4 值 增 幅 极 大 ,确定 土壤 侵蚀 重点 治理 区 域 对 全 面 土壤 侵蚀 防治 具 
关键 词 : 天 山北 坡 ; 水 力 侵蚀 ; CSLE 模 型 ;时 空格 局 ; 


土壤 侵蚀 是 指 在 自然 和 人 为 因素 的 作用 下 , E 
培 及 其 母 质 的 侵蚀 破坏 .运输 和 沉积 的 过 程 。 研 
究 土壤 侵蚀 对 提高 土地 生产 力 , 改 善 环境 、 降 低 风 
沙 ,泥石流 等 自然 灾害 具有 积极 影响 "”。 面 对 日 益 
严重 的 土壤 侵蚀 等 生态 环境 问题 ,迫切 需要 开展 十 
坏 侵 蚀 监 测 工 作 并 探究 驱动 影响 因子 ,对 新 疆 天 山 
北 坡 土壤 侵蚀 综合 治理 提供 数据 支持 ,同时 为 该 区 
域 的 水 土 保持 和 经 济 的 可 持续 发 展 提供 科学 参考 。 

对 于 土壤 侵蚀 的 定量 研究 ,土壤 侵蚀 模型 方法 
是 最 广泛 .最 有 效 的 方法 。 通 用 土壤 流失 方程 (Uni- 
versal Soil Loss Equation, USLE ) 2 和 修正 通用 土壤 
流失 方程 (Revised Universal Soil Loss Equation, 
RUSLE)" 是 国外 经 验 模型 的 典型 代表 。 刘 宝 元 以 
USLE 和 RUSILE 模 型 为 基础 ,经 改进 提出 了 适用 于 
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地 理 探测 器 ; 影响 因素 


全 国 坡 面 侵蚀 的 CSLE 模 型 “。CSLE 模 型 是 目前 发 
展 较为 成 熟 的 土壤 侵蚀 模型 ,其 充分 考虑 到 我 国 陡 
坡 土壤 侵蚀 的 特点 和 水 土 保持 措施 ,其 结构 简单 、 
参数 易于 获取 。CSLE 模 型 已 被 广泛 应 用 于 干 量 区 、 
半 干 旱 区 “探究 土壤 侵蚀 时 空 动 态 变 化 及 影响 
素 。 尽 管 以 上 人 研究 侧重 于 地 理 环境 因素 和 人 类 活 
动 对 土壤 侵蚀 的 影响 ,但 其 并 无 法 量化 影响 因素 对 
土壤 侵蚀 的 影响 程度 ,而 地 理 探测 器 不 仅 可 以 分 析 
单个 自 变 量 对 因 变 量 的 解释 力 ,还 可 以 确定 自 变 量 
之 间 的 交互 作用 对 因 变 量 的 影响 ”。 近 年 来 ,地 理 
探测 带 应 用 广泛 ,其 研究 区 域 大 到 国家 尺度 ,小 到 
乡镇 尺度 ;其 研究 领域 在 社会 科学 、 人 类 健康 自然 
科学 和 环境 科学 等 方面 得 到 应 用 。 部 分 学 者 采用 
地 理 探 测 需 对 土壤 侵蚀 定量 归 因 进行 研究 ”, 可 见 
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将 CSLE 模 型 应 用 于 土壤 侵蚀 定量 评价 并 通过 地 理 
探测 器 探究 土壤 侵蚀 的 空间 分 异 及 其 内 在 驱动 力 ， 
是 区 域 土壤 侵蚀 研究 的 基本 趋势 。 

新 疆 大 部 分 土壤 侵蚀 的 研究 主要 集中 在 降水 
量 相对 较 高 的 伊犁 河谷 区 域 和 ,而 关于 天 山南 北 
坡 的 研究 较 少 , 且 相 关 研 究 主要 集中 在 小 流域 单一 
年 份 土壤 侵蚀 情况 ,缺乏 关于 土壤 侵蚀 的 时 空 动 
态 变化 以 及 定量 归 因 研究 。 天 山 山 区 降雨 存在 明 
显 的 空间 异 质 性 , 北 坡 (迎风 坡 ) 大 于 南 坡 (背风 
坡 )。 天 山北 坡 是 国家 级 水 土 流失 重点 预防 区 , 根 
据 第 一 次 全 国 水 利 普 查 水 土 保 持 情 况 公 报 数据 六 ， 
天 山北 坡 水 蚀 强度 较 大 的 区 域 主要 集中 在 天 山北 
坡 中 段 ,主要 包含 阜 康 市 .乌鲁木齐 市 .昌吉 市 , 呼 
ABEL . 玛 纳 斯 县 石河子 市 . 沙 湾 县 7 个 行政 区 的 
山地 丘陵 区 。 因 此 本 文 以 天 山北 坡 中 段 山 区 为 例 ， 
基于 CSLE 模 型 定量 分 析 2000 一 2018 年 研究 区 的 土 
二 水 力 侵蚀 时 空 变化 特征 ,充分 掌握 土壤 侵蚀 状 
况 ,并 应 用 地 理 探测 器 方法 定量 探究 引起 土壤 侵蚀 
的 多 种 因素 ,为 后 续 水 土 保持 建设 提供 科学 支撑 和 
理论 依据 。 


1 研究 区 概况 与 方法 


1.1 研究 区 概况 

人 研究 区 位 于 天 山北 坡 中 段 ,包含 阜 康 一 沙 湾 共 
7 个 行政 区 所 在 的 山地 区 域 (图 1)。 地 理 位 置 为 
84°50'~88°58'E ,42°55'~44°22'N , 东西 宽 334 km , 南 
北 长 161 km, 总 面积 约 2.80 x 10' km。 地 势 南 高 北 
低 ,海拔 约 为 472~5$251 m。 天 山北 坡 迎 风 背 阳 , 降 
雨量 多 ,太阳 辐射 少 ,较为 湿润 ,多 年 平均 降雨 量 为 
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图 1 研究 区 示意 图 
Fig. 1 Diagram of the study area 


452.03 mm, 降雨 集中 在 5 一 8 月 , 占 年 降雨 量 的 
66.59% 。 年 平均 气温 为 7.2 % ,夏季 平均 气温 为 
23.9% ,冬季 平均 气温 为 -12.5 %。 土 壤 类 型 主要 为 棕 
TEE RSE RG ee Re Mto M 
类 型 多 样 且 时 空 分 布 差异 大 ,主要 由 草本 (72.13% )、 
FEAR (7.10% ) FEAR (0.68% ) 构 成 。 植 被 覆盖 度 在 山 
区 较 高 ,平原 区 较 低 , 在 阴 坡 较 高 , 阳 坡 较 低 。 

1.2 数据 来 源 与 处 理 

依据 新 疆 季 节气 候 及 植被 生长 情况 ,选择 除去 
冬季 外 的 3 一 11 月 为 研究 时 间 段 。 选 取 2000 一 2018 
年 4 一 10 月 的 Landsat 系列 遥 感 影像 数据 ,时 间 分 辨 
率 为 16 d, 空 间 分 辩 率 为 30 m, 数 据 来 源 于 美国 地 
质 调 查 局 USGS (https://www.usgs.gov/) ;从 Earth Data 
(https://earthdata.nasa.gov/ ) 获取 2000—2018 年 MO- 
DIS 数据 产品 MOD13Q1 ,时间 分 辩 率 为 16 d, 空 间 分 
辩 率 为 230 m。 从 中 国 气 象 数据 网 (http:/data.cma. 
cn/) 获 取 2000 一 2018 年 乌鲁木齐 牧 试 站 . 马 鲁 木 齐 
站 及 周边 区 域 共 18 个 气象 站 点 的 日 降水 数据 。 从 
地 理 空间 数据 云 (http://www.gscloud.cn/) 获 取 人 研究 区 
SRTM 空间 分 辨 率 为 30 m 的 DEM 数 据 。 根 据 中 国 
《土地 利用 现状 分 类 》(GB/T 21010-2017)", KAA 
级 分 类 体系 ,结合 研究 区 实际 情况 , 共 划 分 12 类 (天 
然 牧 革 地 、 人 工 牧草 地、 其 他 草地 、 乔 木林 地 、 灌 木 
林地 其 他 林地 .旱地 水 浇 地 、 采 矿 用 地 、 风 景 名 胜 
设施 用 地 .内 陆 滩涂 .农村 宅 基地 ) ,使 用 Landsat 系 
列 遥 感 影像 .高 分 影像 通过 人 机 交互 解 译 获 取 人 研究 
区 2000 年 .2005 年 .2010 年 .2015 年 .2018 年 土地 利 
用 数据 。 本 研究 依据 区 域 土 地 利用 类 型 的 空间 分 
布 特征 ,分别 在 各 期 土地 利用 分 类 图 上 采用 分 层 随 
机 抽样 方法 抽取 200 个 样 点 ,经 野外 调查 和 Google 
Earth 历史 高 空间 分 辨 率 影像 检验 ,总体 分 类 精度 均 
高 于 85% ,满足 研究 所 需要 求 。 从 中 国 科学 院 资源 
环境 科学 与 数据 中 心 (https://www.resdc.cn/) 获 取 研 
究 区 土壤 类 型 数据 。 

MODIS MOD13Q1 数 据 使 用 MRT 软件 进行 格式 
转换 (HDF 转换 成 Tif) ,并 执行 归 一 化 植被 指数 
(Normalized Difference Vegetation Index,NDVI) 数 据 
的 导出 、 投 影 转 换 、 裁 剪 、 最 大 值 合成 法 等 处 理 。 
Landsat 数 据 使 用 ENVI 5.3 软件 进行 预 处 理 ,辐射 定 
标 后 采用 FLAASH (Fast Line-of- sight Atmospheric 
Analysis of Spectral Hypercubes ARAL 进行 大 气 校 
正 。 同 时 ,以 校正 好 的 201$ 年 研究 区 Landsat TM X 
感 影像 为 标准 底 图 ,采用 二 次 多 项 式 拟 合法 对 影像 
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进行 几何 精 纠 正 , SA a BEER NDVI BX HE SF EAT ER 
和 裁剪 等 处 理 。 

野外 调查 采用 分 层 随 机 采样 方法 ,考虑 土壤 类 
型 和 草地 类 型 的 基础 上 ,再 结合 流域 完整 性 及 地 形 
坡 位 .土地 利用 类 型 禁 牧 与 放牧 ,再 根据 样 点 的 均 
匀 性 和 可 达 性 设计 采样 点 。 于 2019 年 7 月 对 82 个 
采样 点 开展 野外 采样 工作 。 样 地 大 小 设置 为 30 m x 
30m ,并 在 对 角 线 的 两 端 及 中 心 位 置 设置 3 个 1 m x 
1 m 的 样 方 。 核 查实 际 土地 利用 类 型 及 采集 0~20 
cm 层 的 土壤 样品 ,采用 重 铬 酸 钾 容 量 法 测定 土壤 的 
有 机 质 含量 ,使 用 激光 粒度 分 析 仪 Microtrac Blue- 
wave S350 测 定 土壤 颗粒 大 小 。 采 用 野外 垂直 向 上 
照相 法 获取 林地 郁 闭 度 及 采用 野外 垂直 向 下 照相 
法 获取 草地 灌木 林地 的 覆盖 度 。 
13 研究 方法 
1.3.1 土壤 水 力 侵 蚀 模 型 ”采用 中 国土 壤 流 失 方程 
(CSLE ) 呈 进行 模型 的 构建 ,计算 公式 如 下 : 

A=RXKXLXSXBXEXT (1) 

式 中 :及 为 降雨 侵蚀 力 因 子 (MJ.mm.hm2 .ha ;天 
为 土壤 可 蚀 性 因子 [thm2.h.hno2.MJ mm; S 分 
别 为 坡 长 因子 坡度 因子 ;BE、7 分 别 为 生物 措施 因 
子 .工程 措施 因子 .耕作 措施 因子 ;4 为 单位 面积 的 
年 平均 土壤 流失 量 (thm?…a')。 将 RKLSBE.T 
因子 的 栅 格 数据 进行 空间 相 乘 , 当 土 地 利用 类 型 为 
耕地 时 ,在 生物 措施 因子 .耕作 措施 因子 两 者 中 , 选 
取 耕 作 措 施 因 子 与 其 他 5 个 因子 栅 格 数据 相 乘 ; 当 
土地 利用 类 型 为 非 耕地 时 , 则 选取 生物 措施 因子 与 
其 他 5 个 因子 栅 格 数据 相 乘 ,得 到 天 山北 坡 中 段 山 
区 土壤 侵蚀 模 数 数据 。 根 据 国 家 水 利 部 颁布 的 《 土 
壤 侵 蚀 分 类 分 级 标准 》(SL190-2007)59 将 土壤 侵蚀 
模 数 分 为 0~200 tkm- a! 200-2500tkm2.a '、 
2500~5000 t-km™?+a™',5000~8000 t- km a"! 8000~ 
15000 t-km™*+a™, > 15000 t-km?+a'FE6 PER 

(1) 降雨 侵蚀 力 因 子 (R) 

通过 对 各 种 算法 的 比较 及 数据 获取 情况 ,采用 
章 文 波 等 "5 根据 日 降雨 量 数据 计算 半月 降雨 侵蚀 
力 的 算法 模型 。 为 了 与 生物 措施 因子 的 时 间 步 长 
保持 一 致 ,将 半月 降雨 侵蚀 力 模型 改进 为 16 4 降雨 
侵蚀 力 模型 。 计 算 2000—2018 年 18 个 气象 站 点 
3—11 月 (以 16d 为 间隔 ) 的 降雨 侵蚀 力 ,并 将 一 年 
中 3 一 11 月 的 降雨 侵蚀 力 累 加 得 到 年 降雨 侵蚀 力 。 
采用 反 距 离 权 重 插 值 (Inverse Distance Weight, 
IDW) 方 法 "“ 计 算得 到 2000 一 2018 年 研究 区 降雨 侵 
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蚀 力 因子 。 

(2) 土壤 可 蚀 性 因子 (K) 

相 较 于 直接 测定 法 和 公式 法 计算 土壤 可 蚀 性 
子玉 值 ,USLE 模 型 中 Wischmeier 等 咏 提 出 的 诺 度 
图 模型 (Nomogragh Model ,Nomo ) 计 算 天 值 , 虽 比较 
繁 天 ,但 精确 度 较 高 ,所 得 的 K 值 较 稳 定 。RUSLE 
模型 中 对 Nomo 模 型 修正 得 到 修正 诺 席 图 模型 
(Modified Nomogragh Model, Mnomo) 。 采 用 Mnomo 
模型 计算 各 土壤 样品 的 K 值 ,得 到 各 土壤 类 型 的 平 
均 K 值 ,并 赋值 于 土壤 类 型 空间 分 布 图 ,得 到 研究 区 
土壤 可 蚀 性 因子 。 

(3) 坡 长 坡度 因子 (L、5) 

采用 Foster 等 所 根据 分 段 坡 的 公式 计算 坡 长 因 
FP) BEE < 10? 选 用 MeCool 等 所 的 公式 ,坡度 > 10° 
选用 Liu 等 2 的 公式 计算 坡度 因子 。 

(4) 生物 措施 因子 (B) 

D ESTARFM 时 空 融 合 模 型 

Landsat 数 据 易 受 云雨 天 气 的 影响 , 且 无 法 获取 
连续 的 高 质量 时 序数 据 ;MODIS 数据 具有 较 好 的 时 
效 性 ,但 空间 分 辨 率 较 低 。 因 此 采用 Zhu 等 ”提出 
的 增强 型 自 适应 反射 率 时 空 融 合 模型 (Enhanced 
Spatial and Temporal Adaptive Reflectance Fusion 
Model, ESTARFM) 用 于 将 Landsat NDVI 空间 数据 
(分 辩 率 30 m, 16 d) 和 MODIS NDVI 时 间 数 据 ( 分 辨 
率 250 m, 16 d) 进 行 时 空 融合 ,时 空 融 合 后 得 到 的 数 
据 与 Landast NDVI 数 据 空 间 分辨 率 一 致 , 与 MODIS 
NDVI 数 据 时 间 序 列 一致 ,统一 使 用 年 积 日 表示 数 
据 的 获取 日 期 , 即 第 65~337 d, 得 到 每 年 各 18 景 空 
间 分 辨 率 为 30 m 时 间 分 辨 率 为 16d4 的 NDVI 数 据 。 

D 植被 覆盖 度 估算 模型 

植被 覆盖 度 的 计算 采用 像 元 二 分 模型 3, 计算 
植被 覆盖 度 的 关键 在 于 NDVI,,,( 纯 植被 端 元 的 ND- 
VI 值 ) 和 NDVI,,( 纯 土壤 端 元 的 NDVI 值 ) 的 确定 ， 
以 往 的 研究 中 将 NDVI, 对 应 于 影像 NDVL、 
NDVI , 对 应 于 影像 NDVI，。 由 于 影像 中 噪声 的 
影响 ,使 得 NDVI 值 过 高 或 过 低 。 为 避免 这 种 误差 ， 
在 确定 NDVI,, 和 NDVI, 时 ,根据 研究 区 主要 植被 
覆盖 类 型 的 NDVI 值 与 对 应 实测 样 点 的 植被 覆盖 度 
变化 曲线 ,把 其 关系 代入 像 元 二 分 模型 中 ,得 到 二 
元 一 次 方程 组 ,进而 计算 得 到 不 同 植被 覆盖 类 型 的 
NDVI,., MI NDVI a 2 。 采 用 时 空 融合 得 到 的 NDVI 
数据 ,以 各 时 期 土地 利用 数据 为 掩 膜 ,将 计算 得 到 
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不 同 植被 覆盖 类 型 的 NDVI.. A NDVI a 代入 像 元 Sa 和 
h™ | 
二 分 模型 中 ,实现 对 植被 覆盖 度 的 估算 ,得 到 每 年 qz1- #1 =1- 3" (2) 
no 


各 18 景 的 植被 覆盖 度数 据 ( 以 16 d 为 间隔 )。 

O 生物 措施 因子 B 值 计算 

采用 化 崇 法 等 的 根 据 生物 措施 因子 与 植被 覆 
盖 度 的 关系 式 ,估算 以 16 d 为 间隔 的 生物 措施 因 
子 。 在 原始 RUSLE 模 型 的 基础 上 ,以 16 d 为 间隔 的 
生物 措施 因子 并 结合 各 时 段 降 雨 侵蚀 力 因子 比例 ， 
计算 得 到 年 内 加 权 平 均 生物 措施 因子 ”|。 

(5) 工程 措施 因子 (E) 

依据 遥感 解 译 和 野外 调查 结果 ,研究 区 内 未 发 
现 4 区 域 水 土 流失 动态 监测 技术 规定 (试行 )》2 中 
水 土 保持 工程 措施 因子 赋值 表 所 列 的 工程 措施 类 
型 。 因 此 ,研究 区 工程 措施 因子 统一 赋值 为 1。 

(6) 耕作 措施 因子 (7) 

根据 《区 域 水 土 流失 动态 监测 技术 规定 ( 试 
行 )》” 中 的 耕作 措施 轮作 措施 赋值 表 进 行 赋 值 ,天 


山北 坡 中 段 山 区 属于 西北 干旱 洪 溉 一 熟 兼 二 熟 区 ， 
北 疆 灌溉 一 熟 填 闲 区 ,按照 各 地 块 进行 赋值 ,耕地 


赋值 为 0.281 ,其 他 赋值 为 1。 

1.3.2 地 理 探 测 器 地 理 探 测 器 是 一 种 统计 学 方 
法 ,包括 因子 探测 器 .交互 作用 探测 器 .生态 探测 器 
和 风险 探测 器 4 个 部 分 ,通过 探测 要 素 的 空间 分 层 
异 质 性 来 揭示 某 原因 背后 的 驱动 力 ”。 因 子 探测 器 
可 用 于 分 析 某 种 驱动 力 是 否 引 起 土壤 侵蚀 空间 格 
局 特征 的 主要 原因 ,分 异性 的 大 小 由 g 值 来 度量 ,可 
识别 出 影响 土壤 侵蚀 的 主导 因素 。 交 互 作用 探测 
器 可 用 于 识别 影响 因子 之 间 的 交互 作用 对 土壤 侵 
刨 空 间 分 布 的 解释 力 。 生 态 探 测 器 可 用 于 比较 两 
个 影响 因子 对 土壤 侵蚀 空间 分 布 的 影响 是 否 存在 
显著 差异 。 风 险 探 测 器 可 用 于 判断 单个 影响 因子 
的 各 层 之 间 土 壤 侵 蚀 量 是 否 存 在 显著 差异 ,并 确定 
土壤 侵蚀 高 风险 区 ”。v 值 计 算 公式 如 下 : 


wy LO z 1.0 

HE 0.8 E 0.8 

w & 

= 0.6 4 四 0.6 
© J=1.41x0.01 

2 0.4 pane $ 0.4 

= 0.2 r=0.90** ga 0.2 

ak e P<0.01 z 

0.0 le. i 0.0 
02 04 06 08 * 0.2 

NDVI 


式 中 ;h=1,…,m 为 变量 或 因子 的 分 层 , 即 分 类 或 分 
K; n, 和 nn 分 别 为 层 h 和 全 区 的 单元 数 ; cy Alo? 分 
INEk MEKENI., SSWA SSTA N 
内 方差 之 和 和 全 区 总 方差 。g 值 区 间 为 [0,1], 其 q 值 
越 大 ,该 影响 因素 对 土壤 侵蚀 空间 分 布 的 解释 力 
越 强 。 

土壤 侵蚀 受降 十 .地 形 土壤、 植被 等 多 种 因素 
的 影响 ,本 文选 择 多 年 平均 降雨 量 , 多 年 平均 植被 
覆盖 度 ,坡度 .土壤 类 型 作为 土壤 侵蚀 影响 因子 ,多 
年 平均 土壤 侵蚀 模 数 作为 因 变 量 , 模 型 运行 之 前 对 
连续 型 变量 进行 分 类 处 理 ”。 


2 结果 与 分 析 


2.1 植被 覆盖 度 估算 模型 构建 及 精度 分 析 

利用 草地 (35 个 ) 灌木 林地 (10 个 ) 实 测 采样 点 
的 植被 覆盖 度 及 乔木 林地 (15 个 ) 实 测 采 样 点 的 郁 
闭 度 与 NDVI 构 建 线性 回归 方程 (图 2)。 草 地 和 灌 
木林 地 的 NDVI 与 植被 履 盖 度 均 呈 极 显 车 正 相 关 
(r=0.90,r=0.96;P<0.01)。 乔 木林 地 的 NDVI 与 植 
被 郁 闭 度 呈 极 显 著 正 相关 (r=0.70,P<0.01)。 通 过 
建立 植被 覆盖 度 实测 值 与 估 测 值 之 间 的 线性 关系 
得 到 散 点 图 (图 3) ,决定 系数 请 为 0.93, 均 方 根 误差 
RMSE 为 0.09。 两 者 的 相关 性 散 点 分 布 相 对 集中 ， 
主要 集中 在 1:1 对 角 线 (y=x) 附 近 , 拟 合 趋 势 线 接近 
于 1:1 对 角 线 ,斜率 K 为 1.04, 总 体 预 测 精度 较 高 。 
2.2 土壤 水 力 侵蚀 时 间 格 局 特征 分 析 

从 表 1 可 以 看 出 ,2000 一 2018 年 天 山北 坡 中 有 段 
山区 土壤 水 力 侵蚀 主要 以 微 度 、 轻 度 侵蚀 为 主 , 分 
别 占 总 面积 的 32.34%~40.87% , 33.36% ~43.01% ;中 


1.0r 
其 (c) 
S 0.8 H 
= 06+ 
y=1.44x+0.06 = y=1.01x+0.14 
R=0.93 K 04 R=0.49 
0.96" S gol r=0.70** 
P<0.01 mot P<0.01 
K 0.0 ——————— 
04 06 08 i 02 04 06 08 
NDVI NDVI 


图 2 NDVI 与 草地 灌木 林地 植被 覆盖 度 及 乔木 林地 郁 闭 度 的 线性 关系 


Fig.2 Linear relationship between NDVI and vegetation coverage of grassland, shrub land and canopy density of arbor forest land 
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y=1.04x-0.05 
R=0.93 
RMSE=0.09 


植被 覆盖 度 估 测 值 


0.2 0.4 0.6 0.8 1.0 
植被 覆盖 度 实测 值 
图 3 植被 覆盖 度 实测 值 与 估 测 值 的 线性 关系 
Fig.3 Linearrelationship between measured and estimated 


vegetation coverage 
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势 ,每 年 分 别 以 26.70 km’, 77.47 km 的 速度 减少 。 
总 体 上 ,土壤 水 力 侵蚀 强度 处 于 增加 。 从 图 4 可 知 ， 
2000—2018 年 天 山北 坡 中 段 山 区 年 均 土壤 侵蚀 模 
数 总 体 呈 上 升 趋势 ,每 10 a 增 加 了 576.93 t-km?-a'. 
其 中 2008 年 的 土壤 侵蚀 模 数 最 低 , 为 1673.49 tekm? 
a ;2016 年 的 土壤 侵蚀 模 数 最 高 ,达到 4496.08 t+ 
km?+a', 2000—2018 年 多 年 平均 土壤 侵蚀 模 数 为 
(2817.63 + 940.76)t: km a's 
2.3 土壤 水 力 侵蚀 空间 格局 特征 分 析 

从 表 2 可 以 看 出 ,2000 一 2015 年 天 山北 坡 中 段 
山区 土壤 侵蚀 面积 呈 增 加 趋势 ,与 2000 年 相 比 土壤 
侵蚀 比重 增加 了 8.2% 。 主 要 是 因为 中 度 .强烈 、 极 


强烈 .剧烈 侵蚀 面积 均 呈 逐渐 增加 趋势 ,土壤 侵蚀 
比重 分 别 从 2000 年 的 7.37% 、3.08% .2.97% 3.02% 
增加 至 2015 年 的 13.84% .7.66% 、6.41% 6.39% , 微 


表 1 2000—2018 年 天 山北 坡 中 段 山区 不 同 土壤 水 力 侵 
蚀 强 度 面积 比例 及 变化 率 
Tab. 1 Proportion and change rate of different soil 
hydraulic erosion intensity area in the middle of the = 土壤 侵蚀 模 数 


northern slope of Tianshan Mountains from 2000 to 2018 WF y= ee ow --- 土壤 侵蚀 模 数 变化 趋势 


面积 百分比 /% 面积 变化 率 / A m 
侵 他 级 别 一 = 4000 
最 大 值 MEL PE 标准 偏差 ma’) E i 
微 度 40.87 32.34 36.99 2.69 -26.70 = 3000 oe 
轻 度 43.01 3336 3837 3.11 -77.47 ee) oe | 
中 度 13.84 7.37 10.40 1.39 22.27 ax 2000 
强烈 766 308 532 117 22.10 
极 强烈 ”729 297 463 1.38 28.84 a 1000 
剧烈 752 157 430 195 30.96 
0 Sa AntNOEOWAHACAANMANOLNM CO 
eooocoocoooco Onset tet east ae a 
ooocooooococjocjeocqctocqco 
NANN N NN N NN N NN N NNN ANAN 


度 次 之 , 占 总 面积 的 7.37%~13.84% ;其 他 侵蚀 强度 
占 比 相对 较 少 ,合计 占 总 面积 的 7.62%~22.47%。 
2000—2018 年 研究 区 各 侵蚀 强度 增 速 最 大 的 为 剧 
烈 侵蚀 ,每 年 以 30.96 km? 的 速度 增加 ;中 度 、 强 烈 、 
极 强 烈 侵 蚀 次 之 ,每 年 分 别 以 22.27 km 22.10 km’, 
28.84 km 的 速度 增加 ; 而 微 度 . 轻 度 侵蚀 呈 下 降 趋 


图 4 2000 一 2018 年 天 山北 坡 中 段 山 区 年 
平均 土壤 侵蚀 模 数 变化 


Fig. 4 Changes of annual average soil erosion modulus in 


the middle of the northern slope of Tianshan Mountains from 
2000 to 2018 


表 2 2000 一 2018 年 天 山北 坡 中 段 山 区 不 同 土壤 水 力 侵蚀 强度 面积 及 比例 
Tab. 2 Area and proportion of different soil hydraulic erosion intensity in the middle of 
the northern slope of Tianshan Mountains from 2000 to 2018 
土壤 侵蚀 
年 份 微 度 轻 度 中 度 强烈 极 强 烈 剧烈 合计 


面积 km 比例 /% 面积 /km” 比 例 /% ”面积 km 比例 /% ”面积 km” 比 例 /% 


面积 km 比例 /% 面积 km 比例 /% ”面积 /km 比例 /% 


2000 11362.58 40.54 
2005 9986.28 35.63 
2010 9747.85 34.78 
2015 9063.53 32.34 
2018 10550.34 37.64 


12054.92 43.01 
11590.47 41.35 
11443.27 40.83 

9350.84 33.36 
10130.22 36.14 


2066.71 7.37 864.03 3.08 831.61 2.97 846.98 3.02 
2882.38 10.28 1393.46 4.97 1184.96 4.23 989.29 3.53 
3278.11 11.70 1683.56 6.01 1085.06 3.87 789.00 2.82 
3879.32 13.84 2145.91 7.66 1795.91 6.41 1791.33 6.39 
2736.21 9.76 1485.57 5.30 1527.08 5.45 1597.43 5.70 


16664.26 59.46 
18040.55 64.37 
18278.99 65.22 
18963.31 67.66 
17476.49 62.36 
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BE . 轻 度 侵蚀 面积 比重 共 减 少 了 17.85%。 而 2015 一 表 3 2000 一 2018 年 天 山北 坡 中 段 山区 不 同行 政 区 的 土 
2018 年 土壤 侵蚀 面积 减 小 ,土壤 侵蚀 状况 略 有 好 壤 侵 蚀 模 数 
ke ,土壤 侵蚀 比重 减少 了 5$.30%。 主 要 由 于 中 度 及 Tab.3 Soil erosion modulus of different administrative 
以 上 侵蚀 面积 减少 ,土壤 侵蚀 比重 从 2015 年 的 regions in the middle of the northern slope of Tianshan 


Mountains from 2000 to 2018 


34.30% I> Æ 2018 4E AY 26.21% , 微 度 . 轻 度 侵 蚀 面 ————————— 
行政 区 面积 km 平均 值 /tkm2.a0) ”标准 偏差 


积 比重 共 增 加 了 8.08%。 由 2000 一 2018 年 天 山北 


乌鲁木齐 市 10207.93 3595.54 1383.32 
坡 中 段 山 区 各 行政 区 的 土壤 侵蚀 模 数 可 知 ( 表 3)， 昌吉 Iasg 3355 i 
近 20 a 年 均 土 壤 侵 蚀 模 数 呈 乌鲁木齐 市 > 昌吉 市 > KRT 2096.09 3256.73 1170.54 
ELRETT > 呼 图 辟 县 > 玛 纳 斯 县 > 沙 湾 县 > 石河 子 呼 图 壁 县 3143.32 2251.31 704.94 
市 。 近 年 来 ,乌鲁木齐 市 旅游 业 的 发 展 不 仅 带 来 了 0 ee 
经 济 增长 ,也 对 旅游 景区 的 生态 环境 造成 了 影响 。 7A oe ee aa 

石河 子 市 68.86 1651.66 648.60 


旅游 活动 对 土壤 、 植 被 的 干扰 ,引起 土壤 可 蚀 性 和 

生物 措施 因子 的 变化 , 易 诱 发 土壤 侵蚀 。 甚 次 为 昌 ”用 ,局 部 地 区 的 土壤 侵蚀 有 所 增加 。 由 2000 一 2018 
吉 市 ,畜牧 业 发 展 较为 旺盛 ; 旱 康 市 次 之 ,是 矿产 资 ”年 天 山北 坡 中 段 山 区 土壤 水 力 侵蚀 强度 空间 分 布 
源 集中 分 布 区 。 随 着 旅游 业 .畜牧 业 和 采矿 业 的 发 ”图 可 知 ( 图 5) ,土壤 水 力 侵蚀 强度 较 大 的 区 域 主要 
展 , 植 被 覆盖 等 阻碍 土壤 侵蚀 发 生 的 因素 减少 作 集中 在 乌鲁木齐 市 达 坂 城区 和 南部 高 山区 。 
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图 5 2000 一 2018 年 天 山北 坡 中 段 山 区 土壤 水 力 侵蚀 强度 分 级 
Fig. 5 Classification of soil hydraulic erosion intensity in the middle of the northern slope of Tianshan Mountains from 2000 to 2018 
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2.4 土壤 水 力 侵蚀 的 影响 因素 分 析 
2.4.1 气候 因素 天 山北 坡 处 于 迎风 坡 ,降雨 量 相 
对 较 大 ,影响 土壤 侵蚀 的 主要 气候 因素 为 降雨 。 由 
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坡度 土壤 侵蚀 模 数 总 数 来 看 ,坡度 > 1$" 的 土壤 侵蚀 
模 数 总 数 为 62.12 X 105ta, 占 总 土壤 侵蚀 比重 的 
78.68%。 可 见 土壤 侵蚀 主要 发 生 在 坡度 大 于 15° 的 


2000 一 2018 年 天 山北 坡 中 有 段 山区 不 同 降 雨量 等 

土壤 侵蚀 模 数 可 知 ( 表 4), 当 降雨 量 在 400~450 mm 
时 ,土壤 侵蚀 模 数 最 大 且 波 动 较 大 ,为 (5159.30 + 
2119.70) t+km?+ a; HOY 250~300 mm 时 ,土壤 侵 


KE, EE EH EBOOK FER i, “EE HT fz 
蚀 性 较 高 ,造成 土壤 侵蚀 强度 较 大 。 

2.4.3 土壤 因素 由 2000 一 2018 年 天 山北 坡 中 段 
山区 不 同 土壤 类 型 的 土壤 侵蚀 模 数 可 知 ( 表 6) ,不 


蚀 模 数 为 (3845.19 + 1378.17)t.km2.a; 当 降 雨量 
在 200~250 mm 时 ,土壤 侵蚀 模 数 最 小 且 波 动 较 小 ， 
J (2193.76 + 740.47 )t-km?+a', 


#4 2000—2018 年 天 山北 坡 中 段 山 区 不 同 降雨 量 的 土 
壤 侵蚀 模 数 
Tab.4 Soil erosion modulus of different rainfalls in the 
middle of the northern slope of Tianshan Mountains from 
2000 to 2018 


降雨 量 /mm 面积 /km” ”平均 值 /t:km*…a') ”标准 偏差 
<200 1658.59 2416.40 882.16 
200~250 13440.73 2193.76 740.47 
250~300 5938.33 3845.19 1378.17 
300~350 3537.69 2712.73 984.20 
350~400 1481.54 2707.38 1037.19 
400~450 1390.05 5159.30 2119.70 
> 450 562.61 3247.77 1435.89 


24.2 地 形 因素 ”坡度 是 地 形 地 貌 的 关键 因素 之 
一 ,也 是 影响 土壤 侵蚀 的 重要 因素 。 由 2000 一 2018 
年 天 山北 坡 中 段 山 区 不 同 坡度 等 级 土壤 侵蚀 模 数 
可 知 ( 表 5), 当 坡度 <5° 时 ,土壤 侵蚀 模 数 最 小 且 波 
动 较 小 ,为 (491.33 + 169.21 )t-km?+a!; 当 坡 度 >35° 
时 ,土壤 侵蚀 模 数 最 大 且 波 动 较 大 ,为 (5881.15 + 
2030.53 )t.km?.a"。 坡 度 对 土壤 侵蚀 模 数 有 正 向 影 
响 , 随 着 坡度 的 增 大 ,土壤 侵蚀 模 数 增 大 。 从 不 同 


#5 2000 一 2018 年 天 山北 坡 中 段 山 区 不 同 坡度 的 土壤 
侵蚀 模 数 

Tab.S Soil erosion modulus of different slopes in the 

middle of the northern slope of Tianshan Mountains from 
2000 to 2018 

WAY 人 十 TEES 
pein) PERE apea m 
<5 6211.81 491.33 169.21 3.05 3.87 
5~15 8143.32 1691.67 561.92 13.78 17.45 


坡度 /(*) 面积 /km? 


15~25 5262.29 3211.45 1078.22 16.90 21.41 
25~35 4512.37 4963.40 1703.94 22.40 28.37 
>35 3879.75 5881.15 2030.53 22.82 28.90 


同 土壤 类 型 抵抗 土壤 侵蚀 的 能 力也 各 不 相同 , 呈 草 
秸 土 > 灰 棕 漠 土 > 棕 钙 土 > RAE > AGED 
fh > 灰 漠 土 > 黑 钉 土 > 灰 褐 土 。 从 不 同 土壤 类 
型 土壤 侵蚀 模 数 总 数 来 看 ,土壤 类 型 为 棕 钙 土 E 
HE . 栗 钙 十 时 ,土壤 侵蚀 模 数 总 数 为 42.02x10' ta， 
占 总 土壤 侵蚀 比重 的 71.65%。 可 见 土壤 侵蚀 主要 
发 生 在 土壤 类 型 为 棕 钙 土 草 千 土 、 栗 钙 土 的 区 域 ， 
主要 是 由 于 标 钙 土 . 栗 钙 土 位 于 海拔 1200 ~ 1600 m 
的 北部 山 前 冲积 忆 及 倾斜 平原 区 域 , 植 被 覆盖 较 
差 ,土壤 可 侵蚀 性 较 高 , 且 土 地 面积 较 大 ; 草 秸 土 位 
于 海拔 3000 mm 以 上 ,其 土壤 侵蚀 模 数 最 大 。 


表 6 2000 一 2018 年 天 山北 坡 中 段 山区 不 同 土壤 类 型 的 
土壤 侵蚀 模 数 
Tab. 6 Soil erosion modulus of different soil types in the 
middle of the northern slope of Tianshan Mountains from 
2000 to 2018 


土壤 面积 平均 值 标准 合计 RAE 

类 型 Ikm? tkm2a) MÆ (10% tea") H% 
Hfi 3547.26 4118.67 1453.32 14.61 24.91 
灰 棕 漠 土 1286.54 2903.56 1247.67 3.74 6.37 
棕 钙 土 。 6526.77 2399.07 841.61 15.66 26.70 
栗 钙 土 。 6204.56 1893.80 657.52 11.75 20.04 
EE 1010.11 1840.50 776.76 1.86 3.17 
BER 3639.99 1635.77 632.97 5.95 10.15 
Ti+. 1515.67 1246.43 417.07 1.89 3.22 
WEE 1961.49 759.18 291.66 1.49 2.54 
灰 褐 土 。 2317.17 733.31 270.80 1.70 2.90 


2.4.4 植被 因素 ”由 2000 一 2018 年 天 山北 坡 中 段 
山区 不 同 植被 覆盖 度 等 级 土壤 侵蚀 模 数 可 知 ( 表 
7) , 随 着 植被 覆盖 度 的 升 高 ,土壤 侵蚀 模 数 降低 。 
当 植被 覆盖 度 < 15% 时 ,土壤 侵蚀 模 数 最 大 且 波 动 
较 大 ,为 (4362.28 + 1680.50)t-km?-a!, HRA 
侵蚀 ; 当 植被 覆盖 度 > 75% 时 ,土壤 侵蚀 模 数 最 小 上 且 
波动 较 小 ,为 (561.44+ 上 290.90)tkm2.a, 说 明 较 高 
的 植被 覆盖 度 能 够 对 土壤 进行 保护 ,防止 发 生 土 壤 
侵蚀 。 从 不 同 植被 覆盖 度 土壤 侵蚀 模 数 总 数 来 看 ， 
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表 7 2000 一 2018 年 天 山北 坡 中 段 山区 不 同 植被 覆盖 
度 的 土壤 侵蚀 模 数 
Tab. 7 Soil erosion modulus of different vegetation 
coverage in the middle of the northern slope of Tianshan 
Mountains from 2000 to 2018 


植被 覆 面积 平均 值 标准 合计 ” 模 数 比 
HJE% /km  ě /(te-km?+a') mÆ (10ta H% 
<15 14147.50 4362.28 1680.50 61.72 77.65 
15~30 3794.06 2071.44 710.12 7.86 9.89 
30~45 2474.79 1596.20 567.95 3.95 4.97 
45~60 3412.61 930.85 376.54 3.18 4.00 
60~75 2744.78 716.53 300.64 1.97 2.47 
>75 1435.80 561.44 290.90 0.81 1.01 


当 植 被 覆盖 度 < 15% 时 ,土壤 侵蚀 模 数 总 数 为 61.72 x 
10° ta, 占 总 土壤 侵蚀 比重 的 77.65%。 可 见 土壤 
侵蚀 主要 发 生 在 植被 覆盖 度 <15% 的 区 域 ,主要 是 
由 于 植被 覆盖 度 <15% 的 土地 面积 较 大 ,为 
14147.50 km ,上 且 土 壤 侵蚀 模 数 较 高 。 
25 土壤 侵蚀 影响 因素 定量 分 析 

由 2000—2018 年 天 山北 坡 中 段 山 区 土壤 侵蚀 
影响 因素 探测 结果 可 知 ( 表 8) ,不 同 影响 因子 对 土 
二 侵蚀 强度 的 解释 力 有 明显 差异 ,g 值 由 大 到 小 依 
次 为 :降雨 (0.49) > 土壤 类 型 (0.17) > 坡度 (0.11) > 
植被 覆盖 度 (0.10)。 其 中 降雨 对 土壤 侵蚀 强度 的 解 
释 力 最 强 , 是 影响 土壤 侵蚀 空间 格局 分 布 的 主要 因 
素 ;其 次 为 土壤 类 型 ,其 土壤 理化 性 质 与 土壤 侵蚀 
密切 相关 。 交 互 作用 探测 器 结果 表明 ,不 同 影响 因 
素 通 过 交互 作用 对 土壤 侵蚀 空间 分 布 的 影响 远大 
于 单个 因素 。 其 中 坡度 门 土壤 类 型 (0.28) 为 非 线 性 
增强 ;其 余 因 子 两 两 之 间 交 互 为 双 因 子 增强 。 降 雨 
与 其 他 因子 之 间 的 交互 作用 明显 大 于 剩余 因子 之 
间 的 交互 作用 。 生 态 探 测 器 结果 表明 ,坡度 和 土壤 
类 型 对 土壤 侵蚀 空间 分 布 的 影响 存在 显著 差异 


表 8 2000 一 2018 年 天 山北 坡 中 段 山区 土壤 水 力 侵蚀 影 
响 因 素 探测 结果 
Tab. 8 Detection results of influencing factors of soil 
hydraulic erosion in the middle of the northern slope of 
Tianshan Mountains from 2000 to 2018 


降雨 坡度 土壤 类 型 ” 植被 覆盖 度 
降雨 0.49 = = - 
坡度 0.54 0.11 - = 
土壤 类 型 0.58 0.28 0.17 - 
植被 覆盖 度 0.55 0.18 0.26 0.10 


(P<0.05)。 通 过 风险 探测 器 并 结合 前 文 的 影响 
素 分 析 识 别 出 高 风险 区 ,结果 表明 土壤 类 型 为 棕 钙 
t 草 秸 土 、 栗 钙 土 ,植被 覆盖 度 小 于 15% 坡度 大 
于 15°* 和 降雨 量 在 400 ~ 450 mm 之 间 的 地 区 为 高 风 
念 侵蚀 区 域 。 


3 讨论 


根据 2013 年 第 一 次 全 国 水 利 普查 水 土 保持 情 
况 公 报 数据 显示 ,新 疆 水 蚀 总 面积 为 8.76 x 10* km’, 
其 中 轻 度 侵蚀 面积 占 比 最 大 ,中 度 侵蚀 次 之 ,分 别 
为 74.06% 和 21.40%; 而 强烈 . 极 强烈 .剧烈 侵蚀 面积 
均 占 比较 少 ,分 别 为 2.92% 1.51% 0.11%"; ÆR 
北 坡 中 段 山 区 微 度 . 轻 度 .中 度 .强烈 . 极 强 烈 .剧烈 
侵蚀 面积 占 比 分 别 为 36.99% 、38.37% , 10.40% 、 
5.32% .4.63% .4.30%; 相 比较 而 言 , 天 山北 坡 中 段 山 
区 侵蚀 强度 略 高 于 新 疆 平 均 水 平 。 众 多 学 者 对 天 
山北 坡 的 研究 较 少 且 研 究 主 要 集中 在 小 流域 单一 
年 份 土壤 侵蚀 情况 ,而 区 域 尺度 长 时 间 序 列 土壤 侵 
蚀 模 数 的 合理 估算 及 空间 分 布 特征 对 于 准确 预测 
土壤 侵蚀 尤为 重要 。 本 研究 充分 考虑 到 降雨 和 植 
被 覆盖 的 季节 动态 性 , 均 以 半月 为 步 长 计算 降雨 侵 
蚀 力 因子 及 生物 措施 因子 。 由 于 生物 措施 因子 的 
大 小 对 降雨 引起 土壤 侵蚀 的 作用 非常 敏感 ,因此 有 
必要 确定 植物 在 不 同 生 长 期 对 降雨 侵蚀 的 影响 。 
Schmidt 等 人 的 结论 表明 ,与 单一 的 年 度 生 物 措施 因 
子 相 比 ,关注 四 季 或 一 年 12 个 月 的 生物 措施 因子 的 
动态 方法 有 助 于 减少 年 度 土壤 侵蚀 计算 中 的 误差 ， 
采用 时 空 变化 的 降雨 侵蚀 力 因 子 和 生物 措施 因子 
的 组 合 将 会 更 动态 更 精准 的 评估 土壤 侵蚀 风险 所。 

利用 地 理 探测 器 对 降雨 地形 土壤、 植被 土壤 
侵蚀 因子 进行 定量 分 析 ,结果 表明 降雨 是 引起 土壤 
侵蚀 的 主要 驱动 因素 ,天 山北 坡 山 区 受 季 节 性 强 降 
水 的 影响 ,水 蚀 较 为 强烈 5。 本 研究 表明 ,坡度 对 土 
壤 侵 蚀 模 数 有 正 向 影响 ,这 与 已 有 研究 规律 一 致 …。 
黑 钉 土 及 灰 褐 土 的 土壤 侵 刨 模 数 最 低 。 黑 钙 土 灰 
褐 圭 分别 是 草原 植被 .森林 灌木 下 的 土壤 ,土壤 较 
为 肥沃 ,植被 覆盖 较 好 , 林 下 草 治世 成 。 林 草 的 覆 
盖 能 够 阻挡 降雨 的 冲刷 ,防止 土壤 侵蚀 的 发 生 ， 
此 林 草 地 土壤 不 易 产 生 土壤 侵蚀 ,质地 较 好 的 土壤 
具有 较 强 的 保水 能 力 , 可 以 促进 植被 生长 ,进而 降 
低 土 壤 可 蚀 性 ,与 本 研究 结果 基本 一 致 中 。 随 着 植 
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被 覆盖 度 升 高 ,土壤 侵蚀 模 数 降低 ,主要 是 因为 良 
好 的 地 表 植 被 能 够 防止 雨滴 飞溅 侵蚀 ,减缓 地 表 径 
流 ,发 达 的 地 下 根系 能 够 回 持 土壤 与 改善 土壤 状 
况 ,防止 径流 侵蚀 。 


4 结论 


本 研究 基于 CSLE 模型 ,应 用 样 地 调查 .地理 信 
息 系 统 与 数理 统计 ,并 结合 地 理 探测 器 方法 ,定量 
分 析 了 2000 一 2018 年 天 山北 坡 中 段 山 区 土壤 水 力 
侵蚀 时 空格 局 特征 ,明晰 降雨 .地 形 .土壤 和 植被 因 
素 对 土壤 水 力 侵蚀 的 影响 。 主 要 得 出 以 下 结论 : 

(1) 时 间 序 列 上 , 2000—2018 年 天 山北 坡 中 段 
山区 土壤 水 力 侵蚀 强度 主要 以 微 度 、 轻 度 侵 刨 为 
主 , 分 别 占 总 面积 的 32.34% ~40.87% 、33.36% ~ 
43.01%; 但 近 20 a RWE . 轻 度 侵蚀 面积 呈 下 降 趋 
势 (-26.70 km -a -77.47 km- a`) ,而 其 他 侵蚀 强 
度 面积 均 呈 上 升 趋势 (22.10~30.96 km? a!) ,土壤 侵 
伺 日 趋 严重 。 空 间 分 布 上 , 2000—2018 年 研究 区 多 
年 平均 土壤 侵蚀 模 数 最 大 的 是 乌鲁木齐 市 ,为 
(3595.54 + 1383.32)t-km?+a!; YW B i rh A 
康 市 ,分别 为 (3321.55 + 1182.19) t kn?» al, 
(3256.73 + 1170.54)tkm a ;而 沙 湾 县 石河子 市 
最 低 , 分 别 为 (1752.93+767.79)t'km?.a (1651.66 + 
648.60)t. km a!o 

(2) 不 同 影响 因素 对 土壤 侵 刨 的 响应 各 不 相 
同 。 降 雨 的 大 小 对 土壤 侵蚀 强度 影响 不 同 , 当 降雨 
量 在 250~300 mm ,400~450 mm 时 ,土壤 侵蚀 模 数 最 
大 ;坡度 对 土壤 侵蚀 强度 影响 较为 明显 , 随 着 坡度 
的 增 大 ,土壤 侵蚀 模 数 增 大 ;不 同 土壤 类 型 对 土壤 
侵蚀 的 抵抗 能 力也 不 同 , 呈 草 秸 土 > 灰 棕 漠 土 > BK 
钙 土 > Ht > eet > BHt> KAE > Be 
土 > 灰 褐 土 ; 较 高 的 植被 覆盖 度 能 够 对 土壤 进行 保 
护 并 防止 发 生 土 壤 侵蚀 。 随 着 植被 覆盖 度 的 增加 ， 
土壤 侵蚀 模 数 降低 。 
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Study on temporal and spatial variation characteristics and influencing 
factors of hydraulic erosion in the middle of the northern slope of 
Tianshan Mountains based on CSLE model 


CHANG Mengdi'*, WANG Xinjun'’, LI Na'’*, YAN Linan'’, MA Ke”, LI Juyan’ 
(1. College of Grassland and Environment Sciences, Xinjiang Agricultural University, Urumqi 830052, Xinjiang, 
China; 2. Xinjiang Key Laboratory of Soil and Plant Ecological Processes, Urumqi 830052, Xinjiang, China; 


3. General Station Soil and Water Conservation and Ecological Environment Monitoring of Xinjiang, Urumqi 


830000, Xinjiang, China) 


Abstract: It is of great significance to ecological early warning and soil erosion prevention on mastering the 
spatial law and driving force of soil water erosion differentiation on the northern slope of Tianshan Mountains, 
which could provide theoretical basis and data support for comprehensive management of regional ecological 
environment on the northern slope of Tianshan Mountains. Taking the mountainous area in the middle of the 
northern slope of Tianshan Mountains as an example, the temporal and spatial pattern characteristics (area, 
intensity and geographical distribution) of soil hydraulic erosion in the study area from 2000 to 2018 were 
quantitatively analyzed by means of field investigation, geographic information system, mathematical statistics 
and geographic detectors, and the internal driving forces of rainfall, topography, soil and vegetation on soil 
hydraulic erosion intensity were explored by means of geographic detectors based on the Chinese soil loss 
equation (CSLE). The results showed that: (1) From 2000 to 2018, the intensity of soil hydraulic erosion in the 
middle of the northern slope of Tianshan Mountains was mainly slight erosion and mild erosion, accounting for 
32.34% 40.87% and 33.36%-43.01% of the total area respectively. In recent 20 years, the area of slight and light 
erosion had a downward trend (—26.70, -77.47 km- a`), while the area of other erosion intensities had an upward 
trend (22.10-30.96 km’: a`’), and the overall soil hydraulic erosion intensity was in an increasing trend. (2) The 
overall soil erosion modulus was Urumqi > Changji City > Fukang City > Hutubi County > Manas County > 
Shawan County > Shihezi City. The spatial distribution of erosion intensity in the middle of the northern slope of 
the Tianshan Mountains was closely related to rainfall, topography, soil, and vegetation. The soil types were 
brown calcareous soil, grass felt soil, and chestnut soil. The area with vegetation coverage less than 15%, slope 
greater than 15° and rainfall in the range of 400-450 mm was a high-risk erosion area. (3) The magnitude of the 
differentiation was measured by the q value in the factor detector. The greater the q value, the stronger the 
explanatory power of the influencing factor to the spatial distribution of soil erosion, the rainfall (0.49) > soil type 
(0.17) > slope (0.11) > vegetation coverage (0.10). Different influencing factors obviously enhanced the spatial 
heterogeneity of soil erosion through interaction, and the coupling effect of vegetation coverage and rainfall 
factors made the q value increase greatly. It is of great significance to determine the key control areas of soil 
erosion for comprehensive soil erosion prevention and control. 

Keywords: northern slope of Tianshan Mountains; hydraulic erosion; CSLE model; temporal and spatial pat- 
tern; geographic detector; influencing factors 


